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^ . 1. Introduction 

The Green function of a quantum Hamiltonian (integral kernel of the resolvent) is one of the 



CN I characteristics whose knowledge usually permits to perform the complete spectral analysis 

^ I and to study various perturbations. The aim of the present communication is to obtain explicit 

I expressions for the Green function for a class of spin-orbit Hamiltonians, namely, for the 

^ ' Rashba and Dresselhaus Hamiltonians with uniform magnetic fields, whose study plays a 
central role in the spintronics [IJ. We use an abstract version of the construction [2J, which 

^ I permits to reduce the problem to the well-known Green functions of the Landau Hamiltonian 

ci I and the Laplacian. 



2. Spin-orbit Hamiltonians 



Below we use the Pauli matrices Ox = \ ^ \^ '^y = \ ■ I ' = I ,^ 1 



1 \ / -i \ / 1 

1 '^>= i '^^= -1 



denote the identity 2 x 2-matrix by Oq. 

We consider Hamiltonians of a charged two-dimensional particle in a uniform magnetic 

field B orthogonal to the plane and take into account the spin-orbit interaction. Let A be 

dA ■ dA 

the magnetic vector potential, i.e. B = -zr^ =— ^. In what follows we use the symmetric 

ox oy 

.By Bx. e 
gauge, A(jc,y) = ~^)- Denote as usual pj :— —ihVj and 11^ := pj — Aj, j = x,y. The 

i Deceased (2 April 2007) 
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Hamiltonian without spin-orbit interaction acts in the spinor space L^(M^,C^) and takes the 

form^o = 1^ — n^o-Q. We are interested in the following two types of spin-orbit Hamiltonians. 

The first one, the Rashba Hamiltonian Hr, is of the form 

Or 

^ ~" ' 2 



\e\h 

where /ig = is the Bohr magneton (mg is the electron mass), is the effective ^- 

InieC 

factor, and (Zr is the real-valued Rashba constant (whose dimension is MLr'T^) expressing 
the strength of the spin-orbit interaction. The second one, the Dresselhaus Hamiltonian, is 
given by 

h ' 2 

and Oq is the real-valued Dresselhaus constant (whose dimension is MLr'T~^) expressing the 
strength of the spin-orbit interaction. 

In what follows we use mostly dimensionless coordinates introduced as follows. Denote 

m^CCj 27thc 
xj := — J = R,D. Furthermore, denote by 4>o the magnetic flux quantum, 4>o := , 

e 

Ik Ik by bx 1 

and let b :— —B and a := — A = (— , — — ). Now by setting k := -p, K := k — a and 
9q 4>o 2 2 n 

introducing the coefficient 7:= we rewrite the above Hamiltonians as Hj = Hj, 

2 lUg 2m* 

J = 0, R, D, where 

Hq = Oq, 

Hr = Ho + IkkUk + ybo^, Ur = o^Ky - OyK^, 
Hj) = Ho + 2xDt/D + ybOz, Ud = OyKy - o^K^. 
In what follows we work with these new normalized Hamiltonians Hj. 



3. Reduction to scalar case 

We start with a simple resolvent identity which is an abstract version of the construction 
from [O and which is of crucial importance in all our considerations. 

For simplicity, for any self-adjoint operator A and a complex number E we use the 
notation i? (A, £■) := (A—E)^^. Now let A be a self- adjoint operator acting inacertain Hilbert 
space, a G M. Denote B ■. = A^ + 2aA. 

LetEe C\spec5, then (B-E)-^ = [{A + - {E + a^)]-^ = [(A + a-7])(A + a + 
T])]^\ where r] = VE + a^ (here and below in this section is a fixed continuous branch 
of the square root on the complex plane C with an appropriate cut). 

If both the numbers r] — a and —t] — a are outside of spec A (in particular, if 3 t] 7^ 0) 
and 7] 7^ 0, then 

{{A + a-ri){A + a + ri)y^ = -3-((A-7] + a)"^ -(A + 77 + a)"^). 
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If, in addition, -77 + a ^ specA and T] + oc ^ specA, then (A - 7] + a)"^ - (A + 7] + a)~^ 
= {A + r\ - a){A^ - [r] - a)^)"^ - (A - 77 - a)(A^ - (r] + a)^)"^. As a result, we arrive at 
the identity 

R{B;E) =l-{{A + ri- a)R{A^; (rj - a)^) - (A - 7] - a)R{A^; (7] + af)) 

= ^ (7?(a2; (77 - -7?(a2; (77 + (1) 
1 

+- (i?(A2; (7] - +i?(A2; (7] + a)^)) . 

Our aim now is to calculate the Green functions for //r and Hj^ using Eq. ([T]). 

Consider an operator V] = Uj + /3jOV, where /3j is a real constant which will be chosen 
later, J = R, D. Using the commutation relation Kj^Ky — KyK^ = ib and the elementary 
properties of the Pauli matrices one easily obtains 



vl = iHo + l5l)oo-ba,, V^ = (Ho + 15^)00 + bo,. 



Therefore, 



H} = + 2>i}V} - /3/c7o, J = R, D, 



for /3r:=^Z>, /3d := ^ 
2xR 2xD 



(2) 
(3) 



3.1. The free case 

We consider here the case without magnetic field, b = 0. Then Eq. ([3]) reads simply as 
Hj = Uj+2h]U] with Uf = Ho = -Aoq. Clearly, spec //j = /(spec t/j) for/(jc) =x^ + 2hjx. 
We see that spect/j = M (see appendix), hence spec//j = [— Xj,+oo) and the spectrum 
contains no eigenvalues. 

Note that the Green function Go{x,y;z) of Hq is known explicitly, 

Go(r,r';z) = i-i^o(v^|r-r'|), 

where ^0 is the McDonald function and ^/x > for x> 0. Hence, by ([T]), the Green function 
Gj, which is the integral kernel of the resolvent (Hj — z)^^, has the form 



G,(r,r';.) = ij 



1 



:(t/j-xj)(ifo(Cj+|r-r'|)-ifo(Cr|r-r'|)) 



+KoiCj+\r-r'\)+Ko{Ci 



± 



(z + x^) ±ixj, 



(4) 



where f/j acts as a differential expression with respect to r. Applying the identity ^o(0 — 
—Ki{t) one concludes that 



Gj(r,r';z) 

with 

G]Hr,r';z)=Gf(r,r';z) 
(^o(C|r- 



G]i(r,r';z) Gj^(r,r';z 



12/ 



G2i(r,r';z) 



Gf(r,r';z) 



1 

4n 



i^-(z + x2) 

r'|)-i^o(Cr|r-r'|))+ifo(C|r-r'|)+ifo(Crlr-r1) 
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for both J = R, D, and 



4'(r,r';z) = J^IA - x ^R+i^i (C+|r-r'|) - i;^K,{i;^\v -v'\) 
4;riJ-(z + x2) |r-r'| 

Gi'(r,r';z) = ^y-y')-'<^-^) _ ^ W^K,[i;^\r -r'\) - C^^i (C^ |r-r'|) 



4;riW-(z + x^) |r-r'| 



and Gj ^ (r, r';z) = G]^ (r', r; z ) . We note that such a representation of the Green function was 
essentially obtained in [|2l for special values of z. 

3.2. The magnetic case 

Using ([3]) one can easily calculate the spectrum of Hj. Namely, spec//j = ^(specVj) with 
g{x) = + 2xj — Pf. The spectrum of Vj can be calculated using the results of appendix. 
Note that the spectrum of Hq consists of the Landau levels, spec//o = {|^l(2« + 1) : n G N}. 
Consider first the Rashba case. By ^ we have spec = { |^| (2n + 1—5' signZ?) +/3^ : nE 



N, 5 = ±1} and, respectively, spec Vr = {±y \b\{2n+ 1 — ^signZ?) +/3^: nGN, ^ = ±1}. 
Hence, the spectrum of //r consists of the Rashba levels, spec//R = {£^{n,s) : n eN, s = 



±1}, £^{n,s) = \b\{2n + l-ssignb)±2>c^^ I5l + \b\{2n + l-ssignb). 

For the Dresselhaus case one has, exactly in the same way, specV^ = {|Z?|(2n+ 1 + 



5signZ7) +/3^ : n G N, 5 = ±1}, specVo = {±y |^|(2n + 1 +*signZ7) +j8^ : n G N, 5 = ±1}, 
and the spectrum of //d consists of the Dresselhaus levels, spec//D = '■ n eN, s = 



±1}, £^{n,s) = \b\{2n+l+ssignb)±2>cjj^l5^ + \b\{2n+l+ssignb). We note that the 
formulas for the eigenvalues were obtained e.g. in [3J by a different method. 

Now let us pass to the calculation of the Green functions. Note that Hq has the following 
Green function: 

Go(r.r';.) = i^r(i-5i^)x 

exp(|(rAr')-^(r-r')Xi-^,.;|(r-r')^ 

where ^ is the confluent hypergeometric function [4J. Clearly, 
.y2 .-1 ( iHo-iz-l5l/^±b))-' \ 



^Ho-{z-l5l/^Tb)) 



where +/— corresponds to R/D. Set now rjj := + Xj + /3j^ and 

Ciib) := (77r±xr)2 + Z.-/32, Cd(^) := (tId ± >^Df - b - [5^. 
By ([T]), we have 

+^ ((n' - (m - + (V/ - (r]j + xj)2)-i) . 
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Passing to the Green function, we obtain 
Gj(r,r,z) = — X 

Go(ry;Cr(&))-Go(r,r';C,+W) 




Go(r,r';Cj-(-Z,))-Go(ry;C+(-&)) 



+ 



(6) 



1/ Go(ry;Cr(&)) + Go(r,r';C/(&)) 
2^ Go(r,r';Cr(-^)) + Go(r,r';C/(-&)) y/' 

Here Uj is considered as a differentiation operator with respect to r. Using the identity 
(TV {a, c,x) /dx = —a^{a + 1 , c + 1 ,x) one can write more expUcit expressions for the Green 
function. Namely, 

Gj(r,r,zj-I 21^^^^,.^^ Gf (r,r';z) 



with 

^11, 



2rij 



Gf(r,r';z) = 



l(Go(r,r';CrW)+Go(r,r';Cj+(^): 

+ (Go(r,r';Cr(-^)) -Go(r,r';C/(-^))) + 



i (Go (r, r'; + Go (r, r'; ^+{-b)) ) 

forbothJ = R,D, 

Gi2(r,r';z) = jZ^I {^{x -x') -i{y -y')) (^^^11^ [Go(r,r'; - Go(r,r'; Cj+(-^))^ 



+ 



Fo(r,r';Cr(-&))-^o(r,r';C/(-^)) 



Gh\ry;z) = \b\ {{y-y)-i{x-x')) 



signZ7+ 1 



Go(r,r';Cr(-^))-Go(r,r';C/(-&)) 



Fo(r,r';Cr(-^))-Fo(r,r';Cj+(-^)) 



andGf (r,r';z) = G]2(r',r;z), where 

^o(ry;z) = ^(^-^)r(i-^) 



^2|^| 
/il; 



X 



4. Renormalized Green functions 

In some applications it is necessary to know the renormalized Green function, namely the 
values Gf"(r,r;z) given by Gf"(r,r;z) = limr,^, [Gj(r,r';z) -5(r,r')l , where S(r,r';z) := 
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— — log |r — r'lcTo is the on-diagonal singularity. Terms of this kind appear e.g. when 
calculating the so-called Wigner 7?-matrix. 

Consider first the case b = 0. We will use the representation (HI) for the Green function. 
Using the expansion 

^oW=l7-n2(y) Mm+l)+log2-logz) 

one easily sees that G]^(r, r';z) are continuous functions vanishing at r = r', which means 
that Gf^ is diagonal. To calculate the diagonal terms we use the equality 



Q{z) := }i^_^^{KoiV^zr) +logr) = i^ivC^) - ^log(-z) +log2). 



Hence, 
Gr{r,r;z) 



2i 
1 

27t 



(e(n-e(r))+^(e(r)+e(n) 



V^(l)-^log(-|)+- 



log 



-(z-fXf) +ZXJ 



2i./-(z + x/; 



(z-fxf) -ixj-l 



Co, 



which is independent of r due to the translational symmetry of the Hamiltonian. 
For the magnetic case (b ^ 0) we use the expansions (H 
(-!)« 



^'{a,n+ l,x 

oo 

+ E 



r(a — n] 



4>(a,n-|- l,x) logx 



{n + l)rr 



- ( V/(a + r) - v/( 1 + r) - v/( 1 + n + r) ) .x;' 



+ ■ 



(n-\)\"^^ (a-n)rx'-" 



r(a) 



r=0 



{l-n)r r\ 



The above expansions clearly show that the off-diagonal terms of Gf^ vanish. To express the 
diagonal terms we use the function 

e(z):=lim(c„(r,r';z)+i^log|r-rl) = -i^(,,(i-^)-2,,(l) + logl^ 



then 



Gr(r,r,z) 



ft^z-xj ( Q{Cf{b))-Q{Ci+{b)) 

1 ( e(cr(^))+e(C/W) 







Q{Cf{-b))-Q{C,+ {-b)) 



+ 











Q{^f{-b))+Q{C+i-b)) 



Appendix. Supersymmetric spectral analysis 



For the sake of completeness, here we are going to prove the following 
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Proposition. Let .W^ he Hilbert spaces, Abe a closed densely defined linear operator 

m A* 
A —m 



from M'l to and m > 0. On © consider the operator L \— \ | . Then 



specL = - Y spec(A4* + m^) U y spec (AM + m^) , (A. 1 ) 

and the same correspondence holds for the eigenvalues. 

Proof. First, it is well known Q that specAA* \ {0} = specA*A \ {0}. Clearly, 



A*A + m2 
AA* + m^ 



Therefore, specL^ \ {m^} = spec (AA* + m^) \ {m^}, and for any A G spec AA* \ {0} = 
specAA* \ {0} at least one of the numbers — vA+^, vX + m^ lies in specL. Let us show 
that actually they both are in the spectrum of L. 

Let A > 0, A G specA*A, then there exists a sequence (^„) with 0„ G domA*A C 
domA such that ||0„|| > 1 and lim(A*A — A)0„ = 0. Denote := A + (a/ATw? — 

Clearly, [ ^" | [ ^ ^ I I ^" K which implies 

V^fill > -^ll^^ll '>'^- (A.3) 
By direct calculation, (L — -\/A + m2)v/„ = (v/A+^-m) ^ ^ . Therefore, 




lim(L— \/A + ) = 0. Together with (IA.3I) this implies vA+m^ G specL. 



To show — vA+w? G specL one has to consider the functions 



A — ( \/A — m) 



A* 
A 




where € domAA* C domA*, ||0„|| > 1, and lim(AA* — A)^„ = and to repeat the above 
steps. To finish the proof of Eq. (lA.ll) it is necessary to study the points ±m. Not that it is 
sufficient to consider the situation when ±m is an isolated point in specL (in particular, an 
eigenvalue of L), as for the continuous spectrum the result follows by taking the closure of 
specL\ {— m, m}. 

For m = 0, Eq. (I A. 21) reads as specL^ = spec AA* U spec A*A, and the conditions G specL 
and G specAA* U specA*A are equivalent. 

Assume m 7^ and m is an eigenvalue of L, then there are G domA and (p G domA* 

with II0II + ||(p|| > 1 such that (L — m)(^| = ( ^ ^ 1=0. Clearly, this implies 

\A^-2mq) I 

A(j)n G domA* and A*A^ = 0. If = 0, then also 9 = 0, which contradicts to the choice. 
Therefore, ^ is an eigenvector of A* A. 
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Assume now that is an eigenvalue of A*A, then there is a non-zero vector e 

i I A0 



domA*A C domA with (A*A0„, (/)„) = \\A(I>\\ = 0. Then (L-m) I ^ ) = f 1=0, from 



which m G specL. 

The relationship between the conditions — m G L and G specAA* can be proved in a 
completely similar way. □ 
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